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CN ' The ladder-rainbow truncation of the set of Dyson— Schwinger equations is used to 

study light mesons. The parameters in the effective interaction are constrained by 
the chiral condensate and /„-; the current quark masses are fitted to m w and ttik- 
The dressed quark propagators are in qualitative agreement with recent lattice- 
QCD results at low q 2 while having the correct perturbative behavior at large 
q 2 . The resulting vector meson masses are within 5% of the experimental values. 

^sj The obtained electromagnetic form factors and strong and electroweak coupling 

constants are also in good agreement with the data. At finite temperature, this 
truncation leads to a mean-field chiral phase transition. The spatial pion mass is 
almost constant below this transition, but rises with T close to and above T c . The 
mass of its chiral partner, an idealized a meson, decreases with T until T c , where 
it becomes degenerate with the pion. 
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1 Introduction 



Our aim is to obtain the hadron mass spectrum and observables such as cou- 
pling constants and form factors from the underlying theory, QCD. The .set 
of Dyson-Schwinger equations [DSEs] form a useful tool for this purpose Era. 
In rainbow-ladder truncation, they have been successfully applied to calculate 
the masses and decay constants of light pseudoscalar and vector mesons Era. 
The dressed-quark propagator, as obtained from its DSE, together with the 
Bethe-Salpeter amplitude [BSA] and the qcp) vertex as obtained from the ho- 
mogeneous and inhomogeneous Bethe-Salpeter equations [BSE] respectively, 
fornjLthe necessary elements for form factor calculations in impulse approxima- 
tionErO. Extensions to other mesons, and calculations of strong and electroweak 
d ' decays and branching ratios are straightforwardLJa. 

At finite temperature chiral symmetry is expected to be restored, and 
quarks (and gluons) become deconfined. Since the properties of the pion are 
closely tied to the dynamical breaking of chiral symmetry, an elucidation of the 
T-dependence of these properties is important; particularly since a prodigious 
number of pions is produced in heavy ion collisions. Encouraged by the success 
of the DSE approach to zero-temperature pion_physics, see in Sec. [2|, we have 
extended this approach to nonzero temperatureu, and this is discussed in Sec. [3| 

The DSE for the renormalized quark propagator in Euclidean space is 






S{ P )- 1 =iZ 2 i> + Z 4 m^) + Z l J 7 -± l g'*D lxv {k)— llx S{q)T a ll {q,p), (1) 



where D /J , v (k) is the dressed-gluon propagator, T"(q;p) the dressed-quark- 
gluon vertex, and k = p — q. The most general solution of Eq. (pi) has the 
form S(p)^ 1 = i^A(p 2 ) + B(p 2 ) and is renormalized at spacelike /? accord- 
ing to A([i 2 ) = 1 and B(p?) = m(/i) with m(p) the current quark mass. The 
notation J represents a translationally invariant regularization of the inte- 
gral with A the regularization mass-scale; at the end of all calculations this 
regularization scale can be removed. 

Mesons are described by solutions of the homogeneous BSE 

r H (p + ,p_) = J **K(p,<r,Q) S(q + )r H (q + ,q_)S(q_), (2) 

where p + = p + r/Q andp_ =p— (l — rf)Q are the outgoing and incoming quark 
momenta respectively, and similarly for q±. The kernel K is the renormalized, 
amputated qq scattering kernel that is irreducible with respect to a pair of qq 
lines. This equation has solutions at discrete values of Q 2 = —m 2 Hl where ran 
is the meson mass. Together with the canonical normalization condition for qq 
bound states, it completely determines Tjj, the bound state BSA. The different 
types of mesons, such as (pseudo-)scalar, (axial-)vector, and tensor mesons, are 
characterized by different Dirac -structures. The most general decomposition 
for pseudoscalar bound states isH 

T PS (k + , k-;P) = 75 [iE(k 2 ; k-P: V ) + PF(k 2 ; k-P; V ) 

+ £ G(k 2 ; k ■ P; r,) + a^ k^P u H(k 2 ; k ■ P- 77)] , (3) 

where the invariant amplitudes E, F, G and H are Lorentz scalar functions of 
k 2 and k ■ P. Note that these functions depend on the momentum partitioning 
parameter r\\ however, physical observables are independent of this parameter. 
The bound state BSA, together with the dressed quark propagato|CS, are 
the necessary elements to calculate the electroweak decay constants Ba and, 
in impulse approximation, strong decays. For describing other electroweak 
processes one also needs the dressed qq^ and qqW vertices. These vertices 
satisfy an inhomogeneous BSE: e.g. the qqj vertex r^(p+,p_) satisfies 

r„(p + ,p_) = Z 2lti + J-^ i K(p,q;Q)S(q + )T fl (q + ,q_)S(q-). (4) 

Solutions of the homogeneous version of Eq. (0) define vector meson bound 
states at timelike photon momenta Q 2 = —m 2 . It follows that T fl (p+,p—) has 
poles at those locations. 



2 Model Calculations 

To solve the BSE, we use a ladder truncation, with an effective quark-antiquark 
interaction that reduces to the perturbative running coupling at large mo- 
mentaerQ. In conjunction with the rainbow truncation for the quark DSE, the 
ladder truncation preserves both the vector Ward-Takahashi identity [WTI] 
for the qq^f vertex and the axial-vector WTI. The latter ensures the existence 
of massless pseudoscalar mesons connected with dynamical chiral symmetry 
breakingH. In combination with impulse approximation, this truncation satis- 
fies electromagnetic current conservation U. 

2.1 Effective interaction 

The ladder truncation of the BSE, Eq. (||), is 

K(p, q; P) -> -G(k 2 ) D f ™(k)%^ ® f 7, , (5) 

where D { ™ e (k = p — q) is the free gluon propagator in Landau gauge. The 
corresponding rainbow truncation of the quark DSE, Eq. (fil), is given by 
K(q,p) -> 7„A72 together with g 2 D^(k) -» G(k 2 )D^ e (k). This truncation 
was found to be particularly suitable for the flavor octet pseudoscalar and vec- 
tor mesons since the next-order contributions in a quark-gluon skeleton graph 
expansion, have a si&nificant amount of cancellation between repulsive and 
attractive correctionsEl 

The model is completely specified once a form is chosen for the "effective 
coupling" Q(k 2 ). We employ the AnsatzQ 



G(k 2 ) _4-K 2 Dk\_ k , /u , , ^ 2 lm T{k 2 ) 



■e- k '/" 



±ln 



(1 + k 2 /A 2 



qcd; 



with 7 m = 12/(33 - 2N f ) and f(s) = (1 - exp j^)/s. The ultraviolet behav- 
ior is chosen to be that of the QCD running coupling a(k 2 ); the ladder-rainbow 
truncation then generates the correct perturbative QCD structure of the DSE- 
BSE system of equations. The first term implements the strong infrared-en- 
hancement in the region_D < k 2 < 1 GeV phenomenologically required til to 
produce a realistic valueEj for the chiral condensate of about (240 GeV) 3 . We 
use m t = 0.5 GeV, r = e 2 - 1, Nf = 4, Aq C d = 0.234 GeV, and a renormal- 
ization scale \i = 19 GeV which is well into the perturbative domain Bu. The 
remaining parameters, u = 0.4 GeV and D = 0.93 GeV 2 along with the quark 
masses, arc fitted to give a good description of the chiral condensate, m^iK 
and fn, see Table [l| Note that all invariant amplitudes of the BSA, see Eq. (||), 



Table 1: Results for the light pseudoscalar mesons in GeV. 
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are needed to satisfy the axial- vector WTI u: IZh is the ratio of the residues 
at the meson pole in the pseudovector and pseudoscalar vertices which appear 
in this WTI, and should be equal to one. From Table [j] it is also clear that 
the pseudovector amplitudes F and G do contribute to observables; further- 
more, ithey dominate the asymptotic behavior of the pion electromagnetic form 
factor EJ. 



2.2 Quark propagator 

Within this model, the quark propagator reduces to the perturbative propaga- 
tor in the ultraviolet region. The mass function M(p 2 ) = B(p 2 )/A(p 2 ) behaves 
for large p 2 like 



M(p 2 ) 



(M€d) 



(7) 



where m q is the renormalization-point-independent explicit chiral-symmetry- 
breaking mass. In the chiral limit the behavior is qualitatively different 



M(p 2 ) 



2ir 2 ln 



(qq)° 



(H*fcD 



.1-7*. 



with (qq}° the renormalization-point-indcpcndcnt chiral condensateLj. At 
loop, it is related to the renormalization-dcpendent condensate 



(<jqY 



-Z 4 N r 



1 d A p 



Tr[5 ch iral(p)] , 



(8) 



one- 



(9) 



via (qqy = (uim/Aqcd) 7 " 1 (qq) - 

In the infrared region both Z(p 2 ) = 1/A(p 2 ) and M(p 2 ) deviate signifi- 
cantly from the perturbative behavior, due to chiral symmetry breaking. Quali- 
tatively, our results are similar to those obtained in recent lattice simulationsli-3, 





Figure 1: The quark functions M{p) (left) and Z(p) (right), compared to lattice dataEJ. The 
lattice data correspond to a bare quark mass of 112 MeV, which is between our up/down and 
strange quark mass. Note that the lattice data for M(p) are unreliable above p ss 1.5 GeV. 



see Fig. (|. The strong enhancement of M(p 2 ) in the infrared region is a char- 
acteristic signal for dynamical chiral symmetry breaking. Also the behavior of 
Z(p 2 ) is characteristic: at large p 2 , Z(p 2 ) — 1 up to logarithmic corrections, 
whereas DSE studies typically find Z(0) to be significantly smaller than oncE3. 
This strongly decreased value of Z(0) is now confirmed by lattice data. The 
main difference between our results and the lattice data is the onset of the 
deviation from the perturbative behavior, both for M(p) and for Z(p); this 
mass scale does depend on details of the effective interaction G(k 2 ). 

2.3 Meson Masses and Decay Constants 

Within the same model, we can now solve the BSE for vector mesons, and cal- 
culate the vector meson masses and decay constants. The p and <f> electroweak 
decay constants are related to the p — 7 and 4> — 7 coupling constants respec- 
tively; the decay constant for the K* can be extracted from T-decay into a 
neutrino and a K* via a M^-boson. The results of our model calculations □ are 
shown in Table 0, and are in reasonable agreement with the data. Corrections 



Table 2: Overview of results for light vector meson masses and decay constants, all in GeV. 
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Figure 2: Pseudoscalar and vector meson masses as function of m q /AQcrj: left, mesons 
with equal-mass constituents; right, mesons consisting of one u/d quark and one quark with 
variable mass m q . The current quark mass m q (u) are evaluated at /i = 19 GeV; realistic u/d 
and s quark masses are indicated by triangles, and the corresponding mesons by plusses. 

from pseudoscalar meson loops are estimated to be smallli3, and will generate 
a nonzero width for the vector mesons. 

Although nature provides us with data for fixed quark masses only, it is 
interesting to see how observables such as meson masses vary with current 
quark mass, in particular for comparison with other nonperturbative methods 
such as lattice QCD. As can be seen from Fig. 0, the pseudoscalar meson 
mass mps behaves like (mps) 2 oc m q in the limit m g — > 0, as expected from 
the Gell-Mann-Oakes-Renner [GMOR] relation. Around the strange quark 
mass, the curve starts to deviate significantly from this GMOR behavior, and 
a more linear dependence of the meson mass on m q begins to emergelia. This 
is in agreement with general expectations that for large m q , meson masses are 
simply proportional to the current quark mass. The vector meson masses show 
a much more linear dependence on m q . Over the current mass range explored, 
our results for both the vector and the pseudoscalar mesons can be fitted by a 
form vim = a + fiyjm q / Aqcd + 7 wI(j/Aqcd, with the parameters of Table [| 



Table 3: Fitted parameters for the mesons mass dependence on m q (in GeV). 
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The decay constants all grow with m q , at least in the range of current quark 
masses explored. In the heavy quark limit, one can show analytically that the 
pseudoscalar and vector decay constants vanish like l/-/m g , as expected from 



heavy quark effective theoryQ. The turn-over point for this behavior is larger 
than three times the strange quark mass, at least in the present model. The 
details of this turn-over do of course depend on the details of the effective 
interaction. However, the qualitative behavior, both of the meson masses and 
of the decay constants, is expected to be model-independent. 



2.4 Meson Form Factors 

Meson electromagnetic form factors in impulse approximation are described 
by two diagrams, with the photon coupled to the quark and to the antiquark 
respectively. We can define a form factor for each of these diagramsO, e.g. 



2P v F ail {Q 2 )=N c 



d 4 q 



Tr[S a (q)rZl(q,q + ;P_) 



xS b (q + )iTl(q + ,q_)S b (q_)f;l(q_,q;-P + )] , (10) 



where q — k + \P ■, q± 



hP±hQ, p± 



P ± -iQ. We work in the isospin 



3 ± uss 

and 



symmetry limit, so for the pion form factor we have F 7r (Q 2 ) — F U u U {Q 2 )- The 
charged and neutral kaon form factors are given by F^+ — %F uSu 
and F K o — —-^Fdsd + hFdss respectively. Our results for Q 2 F 7T , Q 2 F K + 
Q 2 F K o are shown in Fig. 13, with the charge radii are given in Table 0. 

Up to about Q 2 = 3 GeV 2 , our results for F v and Fk can be fitted quite 
well by a monopole. Asymptotically, these functions behave like Q 2 F(Q 2 ) — > c 
up to logarithmic corrections. However, numerical limitations prevent us from 



3 GeV 



our result for 
2 



accurately determining these constants.. Around Q 2 

Q 2 F 7r is well above the pQCD result El 16Trf 2 a s (Q 2 ) - 0.2 GeV^, and clearly 
not yet asymptotic. If we continue the calculation for the electromagnetic form 
factor into the timelikc region, we find a pole at the mass of the vector meson 
bound states. Using the behavior of F U u U and F uSS around this pole, we can 



Table 4: Overview of calculated meson charge radii in fm , with o*»t. 
coupling constants, with expt. values extracted from decay widthsl 16 !. 
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Figure 3: The electromagnetic form factors of tho4ight pseudoscalar mesons[J The pion data 
shown here are from an analysis by Bebek et r&k E3- Our results are in good agreement with 
preliminary new data from JLab for the pion^j; new results for F K + are also expected PH. 



extract the coupling constants g p7T n and g<j,KK respectively c3, which govern 
the strong decays p — ► hit and cj) — > KK . The results from this analysis are 
also given in Table [|, and are reasonably close to the experimental data. 

The impulse approximation for the 7*^7 vertex with 7* momentum Q is 






(11) 



N c f d 4 q 



w- 



Tr [S(q) iT v (q, q') S(q') iT^q', q") S{q") T^q" , q; P)] 



where the momenta follow from momentum conservation. In the chiral limit, 
the value at Q 2 — 0, corresponding to the decay 7r° — > 77, is given by the axial 
anomaly and its value g° 77 = \ is a direct consequence of only gauge iavari- 
ance and chiral symmetry; this value is reproduced by our calculations L3 and 
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Figure 4: The 7* 717 form factorE3, with data from CLEO and CELLOE]. 

corresponds well with the experimental width of 7.7 eV. In Fig. |4| we show our 
results with realistic quark masses, normalized to the experimental i g W77 . The 
large-Q 2 behavior of i 7 ' 7 * Tr7 lies in between monopoles fitted to the Brodsky- 
Lepage pQCD asymptotic limit E3, 871- 2 / 2 , and the limit obtained fro m a naive 
analysis of the impulse approximation within the DSE approach EEro, ^7r 2 / 2 . 
Our results for i 7 ' 7 * 7r7 aie compatible with an estimate for the asymptotic be- 
havior by Anikin et al. Ej. 

3 Finite Temperature Extension 

In the Matsubara formalism, the quark propagator can be written as 
S~ x {p,w k ) = n- pA(p 2 ,u 2 k ) +i- fi u k C{p 1 ,ul) + B(p 2 ,ul) 



Z 2 if • p + Z 2 174 uj k + Zim{n) + S(p, Wfe) . 



(12) 
(13) 



where w k = (2k + 1) nT is the fermion Matsubara frequency and T the tem- 
perature. The functions A and C are renormalized such that A(p 2 ,o> 2 ) = 1 = 
C(p 2 ,uj 2 ), and B(p 2 ,ojq) = m(fi) at p 2 + uj 2 = fi 2 . The approximations we 
use are similar to the T = studies discussed in the previous sections. The 
rainbow-ladder truncation for the self-energy is 

' A d 3 c 



£(ftwfc) 



r£ 



i=-c 



(2tt) 



:9 D^if- q, fi fe _ i )7 M S'(q, w,^ , 



(14) 



where g 2 D tlUl the effective interaction or effective gluAO-nxopagator, at finite T 
is taken to be the minimal extension of a T = modelBEillEl with all parameters 
fixed at T = 0. This truncation leads to a second-order chiral phase transition 

(qq)°<x(l-T/T c f, T/T C U, (15) 

where (3 is the critical exponent for chiral symmetry restoration. Independent 
of the details of the effective interaction, we find mean field behavior, j3 = 
1/2, in rainbow-ladder modelsE3. However, one has to go to extremely small 
current quark masses, several orders of magnitude smaller than realistic u/d 
quark masses, to observe this mean-field behavior. The critical temperature 
is sensitive to details of the model, but with the model parameters fitted to 
T = pion observables, one typically finds T c ~ 0.15 GeV. 

3. 1 Meson correlations at finite temperature 

In order to study meson correlations u we solve the inhomogeneous scalar and 
pseudoscalar vertices T S ^ PS for the lowest external (meson) Matsubara mode, 
Sl = 0. Bound states will appear as poles in these vertices; for this choice of 
the external momentum these pole&xprrespond to the spatial modes, and thus 
the masses to the spatial masses oH, which are the type of masses that are 
usually calculated in lattice simulations. 

For to = the T = inhomogeneous BSEs exhibit poles at m v = and 
to ct = 0.56 GeV, with m a = 0.59 GeV at m = 5.7 MeV. This low-mass scalar is 
typical of the rainbow-ladder truncation, although there is some model sensi- 
tivity. The rainbow-ladder truncation yields degenerate isoscalar and isovector 
bound states, and ideal flavor mixing in the 3-flavor case; improvements be- 
yond the ladder truncation are required in order to describe the observed scalar 
mesons. Furthermore, in the isoscalar-scalar channel, it will be necessary to 
include couplings to the dominant nn modeo because of the large a width. In 
the absence of such corrections, the a properties elucidated herein are strictly 
only those of an idealized chiral partner of the it. 

These poles evolve with T in such a wayl3 that at the critical temperature, 
we have degenerate, massless pseudoscalar and scalar bound states, see Fig. g| 
Above T c , these bound states persist, becoming increasingly massive with in- 
creasing T and approaching 2nT, the thermal mass of two free quarks. Similar 
features are also observed in numerical simulations of lattice-QCDE3. Above 
T c all but the leading Dirac amplitude vanishes and the surviving pseudoscalar 
amplitude is pointwise identical to the surviving scalar one. These results in- 
dicate that the chiral partners are locally identical above T Cl they do not just 
have the same mass. Approaching T c from below, our analysis yields for the 
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Figure 5: The 7r and cr mass, decay constant ifis function of T, both in the chiral limit (left) 
and with realistic current quark mass (right) tj. 



a mass m a oc (1 — T/T c )@ within numerical errors. Note that this is the same 
behavior as the chiral condensate and other, equivalent order parameters for 
chiral symmetry restoration, such as the decay constant /^. 

For nonzero current-quark masses chiral symmetry restoration is exhibited 
as a crossover rather than a phase transition. That the transition has become 
a crossover is evident in the behavior of f n , see the right panel of Fig. |5[ The a 
mass exhibits a dip around the phase transition. The scalar and pseudoscalar 
meson masses become indistinguishable at T ~ 1.2 T c . 

With the BSAs obtained from the homogeneous BSE, we can now calculate 
the T-dependence of other meson properties u. The T-dependence of 17^77 is 
calculated in a similar fashion as at T = in the previous section. At T — this 
coupling saturates the Abelian anomaly: 311-77 = h-. The interesting quantity is: 
g® 7//^) which vanishes at T c in the chiral limit Ej: g° 77 goes to zero linearly. 
Thus, in the chiral limit, the coupling to the dominant decay channel closes 
for both charged and neutral pions. For m ^ both the coupling: 3^77/ S-k-, 
and the corresponding width exhibit a crossover with a slight enhancement in 
the width as T — > T c due to the increase in m r . 

Also the isoscalar-scalar-7T7r strong coupling vanishes at T c in the chiral 
limit, which can be traced to i? c hirai —> 0- F° r r ^ 7^ 0, the coupling reflects the 
crossover. However, that is irrelevant because the width vanishes just below T c 
where the a meson mass falls below 2m„- and the phase space factor vanishes. 
A similar phenomenon occurs for the pirir decay width: with a separable model 
for the effective interaction it was shown i3 that slightly above T c the spatial 
mass of the p becomes smaller than 2tti w , and thus the strong decay of a p into 
pions becomes impossible because of phase space consideration. However, one 
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should keep in mind that these are the spatial masses, and that effects such as 
collisional broadening are not yet taken into account. 
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